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ABSTRACT 



>■ ' Received: 1 1 August 2006 / Accepted: 27 March 2007 

i> : 

G\ . 
(N ■ 

I . Aims. We observed the high redshift (z = 3.969) GRB 050730 with Swift and XMM-Newton to study its prompt and afterglow emission. 
' Methods. We carried out a detailed spectral and temporal analysis of Swift and XMM-Newton observations. 
CD Results. The X-ray afterglow of GRB 050730 was found to decline with time with superimposed intense flaring activity that extended 
over more than two orders of magnitude in time. Seven distinct re-brightening events starting from 236 s up to 41.2 ks after the burst were 
' observed. The underlying decay of the afterglow was well described by a double broken power-law model with breaks at ti = 237 ± 20 s and 
U = 10.1 ^2 2 ks. The temporal decay slopes before, between and after these breaks were a\ = 2.1 + 0.3, a 2 = 0.44*^ and = 2.40^*, 
[ ■ respectively. The spectrum of the X-ray afterglow was well described by a photoelectrically absorbed power-law with an absorbing column 
K> ' density A^=(l. 28^25) x 10 22 cirT 2 in the host galaxy. Evidence of flaring activity in the early UVOT optical afterglow, simultaneous with that 
^ observed in the X-ray band, was found. Strong X-ray spectral evolution during the flaring activity was present. The rise and decay power-law 
' slopes of the first three flares were in the range 0.8-1.8 using as zero times the beginning and the peak of the flares, respectively. In the majority 
of the flares (6/7) the ratio Ar/f p between the duration of the event and the time when the flare peaks was nearly constant and ~ 0.6-0.7. We 
showed that the observed spectral and temporal properties of the first three flares are consistent with being due both to high-latitude emission, 
as expected if the flares were produced by late internal shocks, or to refreshed shocks, i.e. late time energy injections into the main afterglow 
shock by slow moving shells ejected from the central engine during the prompt phase. The event fully satisfies the E p -Ei io Amati relation 
while is not consistent with the E p vs. £j el Ghirlanda relation. 
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1. Introduction 

The successful launch on 2004 November 20 of the Swift 
Gamma-ray Burst Explorer dGehrels et al. 2004b has opened 
a new era in the study of Gamma Ray Bursts (GRBs). The 
autonomous and rapid slewing capabilities of Swift allow the 
prompt (1-2 minutes) observation of GRBs, discovered and 
localised by the wide-field gamma-ray (15-350 keV) Burst 



Alert Telescope (BAT, Barthelmy et al. 2005a 1, with the two 



co-aligned narrow-field instruments on-board the observatory: 
the X-Ray Telescope (XRT, IBurrows et al. 2003al) . operating 
in the 0.2-10 keV energy band, and the Ultraviolet/Optical 
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Telescope (UVOT, Roming et al. 2005), sensitive in the 1700- 
6000 A band. The Swift unique fast pointing capability is cru- 
cial in the X-ray energy band, where the reaction times of other 
satellites are limited to time scales of several hours. With Swift, 
thanks also to the XRT high sensitivity, it is possible for the first 
time to study in detail the evolution of the X-ray afterglows of 
GRBs during their very early phases. 
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Indeed, one of the main results of Swift is the identification 
of unexpected and complex features in the early X-ray 
afterglows. In particular, three distinct phases are observed in 
the majority of GRB light curves: an early (f<500 seconds 
from the trigger) steep decline, with a power-law index of ~ 3, 
a second (f<10 ks) very shallow phase with a slope of ~ 0.5, a 
third phase characterized by a more conventional decay slope 
of ~ 1 (e.g. Tagliafer ri et al. 2005| ICusu mano et al. 2006a, 
INousek et al. 2006110' Brien et al. 20061 In a few cases (GRB 
050315, |Vaughan et aL2006l GRB 050 318, IPerri et al. 20031 
GRB 050505, |Hurke tt et al. 20061 GRB 05052 5A, 
IBlustinetal. 20061 GRB 060614, |Mangano et al. 2007| > a 
further steepening with a decay slope of ~ 2, consistent 
with a jet break, is observed. Moreover, Swift had detected 
in about one half of the bursts strong flaring activity in 
the X-ray energy band superimposed on the afterglow 
decay (e.g. IBurrows eT al. 2005b; Romano et al. 2006, 
IFalcone et al. 20061 

The understanding of the origin of these bright X-ray 
flares is intensively discussed in the literature. A mech- 
anism proposed as responsible for the flaring activity is 
late internal shocks (IBurrows et al. 2005bl IFan & Wei 2 005 ; 



Zhang et al. 2006 Liang et al. 20061. In this scenario the X 



ray flares are produced by the same internal dissipation pro- 
cesses which cause the prompt emission, likely internal shocks 
within the expanding fireball occurring before it is decel- 
erated in the external medium (e.g. Rees & Meszaros 1994). 
This model requires that the GRB central engine is still ac- 
tive after the end of the prompt emission and various mech- 
anisms providing such extended internal activity have been 
put forward (e.g. |King et al. 2"005] IPerna et al. 2006b . An al- 
ternative scenario has been recently considered by Guetta et 
al. (2007) who, based on a detailed analysis of the X-ray 
flaring activity observed in the afterglow of GRB 0507 13 A, 
interpreted the X-ray flares as due to refreshed shocks, i.e. 
late time collisions with the main afterglow shock of slow- 
moving shells ejected from the central engine during the 
prompt phase {Regs & Mes zaros 19981 IKumar & Piran 20001 
ISari & Meszaros 20001) . 

In this paper we present a detailed analysis of Swift and 
XMM-Newton observations of GRB 050730, focusing on the 
intense and extended X-ray flaring activity that characterizes 
its afterglow. In Section [2] the observations and the data re- 
duction are presented, in Section [3] we describe the tempo- 
ral analysis and Section H] is dedicated to the spectral anal- 
ysis. The results are discussed in Section [5] and finally in 
Section [6] we summarize our findings. Throughout this pa- 
per errors are quoted at the 90% confidence level for one pa- 
rameter of interest (Ax 2 = 2.71) unless otherwise specified. 
We adopted the standard ACDM cosmological parameters of 
Q m = 0.27, Q A = 0.73 and H = 70 km s" 1 Mpc" 1 . Times are 
referenced from the BAT trigger To while temporal and spec- 
tral indices are written following the notation F{t, v) oc t~ a v~^. 
Results on the optical spectrum of the afterglow of this 
GRB are reported by|Starling et al. (20"05)l|Chen et al. (2005a)| 




100 
Time since BAT trigger (s) 



and |Prochaska et al. (2006)[ Multi- wavelength observations of 
the afterglow of GRB 050730 are presented by Pandey et 
al. (2006). 



Fig. 1. BAT 20-150 keV background subtracted light curve of 
the prompt emission of GRB 050730. Data are binned to 10 
seconds resolution and errors are at the lcr level. The horizontal 
dashed line indicates the level. 



2. Observations and data reduction 

2.1. Swift BAT 

The BAT detected and located GRB 050730 at 
r =19:58:23 UT on 2005 July 30 (Holland et al. 20051 ). 
On the basis of the refined ground analy- 
sis dMarkwardt et al. 2005]) . the BAT position is 
RA(J2000)=212?063, Dec(J2000)=-3°740, with a 90% 
containment radius of 3'. The BAT prompt emission light 
curve (Fig.[TJ is characterized by a duration of T90 = 155+20 s. 
The emission starts ~60 s before the trigger, peaks at ~ 10 s 
after the trigger and declines out to ~ 180 s after the trigger. 

The time averaged (over T90) spectrum in the 15-150 keV 
energy band is well described by a power-law model with en- 
ergy index /3baj = 0.5 + 0.1 and^ = 0.71 (with 56 degrees 
of freedom, d.o.f.). The total fluence in the 15-150 keV band 
is (2.4 ± 0.3) x 10~ 6 erg cirT 2 . Assuming as redshift z = 3.969 
(IChen e t al. 2005b , see Sect. 12. 21 , the isotropic -equivalent radi- 
ated energy in the BAT bandpass (74.5-745.4 keV in the burst 
rest frame) is E®f = (8.0 ± 1.0) x 10 52 erg. 

2.2. Ground-based Observatories 



Following the identification of the optical counterpart by 
UVOT (Holland et al. 20051 ). the field of GRB 050730 was ob- 
served by numerous ground-based telescopes. The afterglow 
detections in the R dSota et al. 20051 ). r' (IGomboc et al. 20051 ). 
/ and J (Cob b et al. 20051 ) bands were soon distributed via the 
GRB Circular Network (GCN). Observations in the optical 
band with the MIKE echelle spectrograph on Magellan II led to 
the GRB redshift measurement z = 3.969 (IChen et al. 20 05b) 
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based on the detection of a strong hydrogen Lya absorption line 
and of several absorption lines from other elements. The red- 
shift was later confirmed with observations by the ISIS spectro- 
graph on the William Herschel Telescope dRol et al. 2 005 ), the 
IMACS imaging spectrograph on the Magellan Observatory 
Baade Telescope dHolman et al. 2005al . FORS1 and UVES on 
the Very Large Telescope dD'Elia et al. 2 005 ). 

The optical afterglow decay at later times was followed 
by several optical telescopes. In particular, measurements 
were made in the R band up to ~4 days after the trig- 
ger dHolman et al. 2005al 2005b; IBurenin et al.~2005l 
IKlotz et al. 20531 |Damerdji et al. 2005| ID'Elia et al. 20051 



IBhatt & Sahu 20051 [Kannapp an et al. 2005) . 

Finally, a radio afterglow with a flux density 
F T = 145 + 28 //Jy at 8.5 GHz was detected about 2 days 
after the trigger using the Very Large Array (Cameron 2005). 

2.3. Swift UVOT 



Swift UVOT began to observe the field of GRB 050730 
at 20:00:22 UT, 119 seconds after the BAT trigger. The first 
100 seconds exposure in the V band led to the identifi- 
cation of the optical afterglow at RA(J2000)=14 h 08 m 17 s .09, 
Dec(J2000)=-03°46'18'/9 dHolland et al. 20051 1. 

We refined the preliminary photometric analysis 
dBlustin et al. 20051 1 processing all UVOT data using the 
standard UVOT software package (Swift software v. 2.1 in- 
cluded in the HEAsoft package v. 6.0.2). The flux in all filters 
was estimated by integrating over a 3.5" region. A background 
region for subtraction of the sky contribution to the flux has 
been selected in a relatively empty part of the field of view. 
The results are listed in Table [1] No significant detection was 
found in the U and UV filters, which is consistent with the 
high redshift measured for this GRB. All the magnitudes 
are corrected for Galactic extinction (E(B - V) = 0.049, 
Schlegel et al."T998) . 



2.4. Swift XRT 

The XRT observations of the GRB 050730 field started at 
20:00:28 UT, 125 seconds after the BAT trigger, with the in- 
strument in Auto State. After a first exposure in Image Mode 
(see IHill et al. 20041 for a description of readout modes) dur- 
ing which no on-board centroid was determined, the instrument 
switched into Windowed Timing (WT) mode for the entire first 
Swift orbit from 7o+133 s to Tq+794 s. Starting from the sec- 
ond orbit (7o+4001 s), the instrument was in Photon Counting 
(PC) mode for 26 consecutive orbits until 1 1 :49:22 UT on 2005 
August 1 (To+143.4 ks). The field of GRB 050730 was re- 
observed with the XRT from August 3 starting at 15:28: 1 1 UT 
(r +329.4 ks) until August 5 13:57:02 UT (J +496.7 ks). 

The XRT data were processed with the XRTDAS soft- 
ware (v. 1.7.1) included in the HEAsoft package (v. 6.0.4). 
Event files were calibrated and cleaned with standard filter- 
ing criteria with the xrtpipeline task using the latest calibration 
files available in the Swift CALDB distributed by HEASARC. 



Events in the energy range 0.3-10 keV with grades 0-12 (PC 
mode) and 0-2 (WT mode) were used in the analysis (see 
IBurrows et al. 2003al for a definition of XRT event grades). 
After the screening, the total exposure time for the first XRT 
observation was 649 seconds (WT) and 58480 seconds (PC), 
while for the follow-up observation the PC exposure time was 
34669 seconds. 

In the 0.3-10 keV PC image of the field a previ- 
ously uncatalogued X-ray source was visible within the 
BAT error circle with coordinates RA(J2000)=14 h 08 m 17.2, 
Dec(J2000)=-03°46'19". This position, derived using data 
not affected by pile-up (orbits 5-26, see below), has 
a 90% uncertainty of 3.5" using the latest XRT bore- 
sight correction (Mor etti et al. 20 06) and is consistent with 
the position of the optical counterpart (Holl and et al. 2 005; 



Jacques & Pimentel 2005 1. 



For the WT mode data, events for temporal and spectral 
analysis were selected using a 40 pixel wide (1 pixel corre- 
sponds to 2.36") rectangular region centered on the afterglow 
and aligned along the WT one dimensional stream in sky co- 
ordinates. Background events were extracted from a nearby 
source-free rectangular region of 50 pixel width. For PC mode 
data, the source count rate during orbits 2-A was above ~ 0.5 
counts s _1 and data were significantly affected by pile-up in the 
inner part of the Point Spread Function (PSF). After comparing 
the observed PSF profile with the analytical model derived by 
Mor etti et al. (2005)] we removed pile-up effects by excluding 
events within a 5 pixel radius circle centered on the afterglow 
position and used an outer radius of 30 pixels. From orbit 5 
the afterglow brightness was below the pile-up limit and events 
were extracted using a 10 pixel radius circle, which encloses 
about 80% of the PSF at 1.5 keV, to maximize the signal to 
noise ratio. The background for PC mode was estimated from 
a nearby source-free circular region of 50 pixel radius. Source 
count rates for temporal analysis were corrected for the fraction 
of PSF falling outside the event extraction regions. Moreover, 
the loss of effective area due to the presence of 2 CCD hot 
columns within the extraction regions was properly taken into 
account. The count rates were then converted into unabsorbed 
0.3-10 keV fluxes using the conversion factor derived from the 
spectral analysis (see Sect. [4]). 

For the spectral analysis, ancillary response files were gen- 
erated with the xrtmkarf task applying corrections for the PSF 
losses and CCD defects. The latest response matrices (v. 008) 
available in the Swift CALDB were used and source spectra 
were binned to ensure a minimum of 20 counts per bin in order 
to utilize the^ 2 minimization fitting technique. 



2.5. XMM-Newton 

XMM-Newton follow-up observations of GRB 050730 
started 26.4 ks (for the two EPIC-MOS cameras) and 29.4 ks 
(for the EPIC-PN) after the initial BAT Trigger. The XMM- 
Newton ODF (Observation Data Files) data were pro- 
cessed with the epproc and emproc pipeline scripts, using 
the XMM-Newton SAS analysis package (v. 6.5). A bright 
rapidly decaying source was detected near the aim-point of 
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Fig. 2. Swift XRT (filled circles) and XMM-Newton (filled squares) 0.3-10 keV light curve of the X-ray afterglow of GRB 050730. 
The BAT 20-150 keV prompt emission light curve, extrapolated to the 0.3-10 keV band, is also shown as filled triangles. UVOT 
optical data in the V and B bands, arbitrarily scaled for comparison with the X-ray band, are indicated with open squares and open 
triangles, respectively. The solid line is the best fit model to the XRT and XMM-Newton light curve. The dashed line represent 
the underlying double broken power-law decay (see Sect. 13. II ). The dotted line is a power-law model with temporal decay index 
ay = 0.3 normalized to fit data points in the V band. Errors are at the lcr level. 



all three EPIC detectors and the afterglow was localized at 
RA(J2000)=14 h 08 m 17 s .3, Dec(J2000)=-03°46'18'/5. The dura- 
tion of the XMM-Newton follow-up observation was 33.7 ks 
(MOS) and 30.4 ks (PN). After screening out times with high 
background flaring, the dead-time corrected net exposures were 
25.0ks (MOS) and 17.9 ks (PN). All three EPIC cameras (PN 
and 2 MOS) were used in Full Window Mode, with PN and 
MOS2 using the "Thin" filter and MOS1 using the "Medium" 
optical blocking filter. 

Source spectra and light curves for all 3 EPIC cameras were 
extracted from circular regions of 30" radius centered on the 
afterglow. Background data were taken from a 60" radius cir- 
cle on the same chip as the afterglow, but free of any back- 
ground X-ray sources. The data were further screened by in- 
cluding only good X-ray events (using the selection expres- 
sion flag=0 in evselect), by including events with single and 
double pixel events (pattern<=4) for the PN and by selecting 
single to quadruple pixel events (pattern<=12) for the MOS. 
Data below 300 eV and above lOkeV were also removed. 

For the temporal analysis, we adopted the light curve from 
the MOS data, primarily because the PN data was heavily af- 



fected by background flares towards the end of the observation, 
while the MOS covered a wider duration at the beginning of the 
observation. The data from the two MOS detectors were com- 
bined and the count rate to 0.3-10 keV unabsorbed flux conver- 
sion factor was calculated from the best fit absorbed power-law 
spectrum (see Sect.Hl). 

For the spectral analysis, ancillary and redistribution re- 
sponse files for fitting were generated with the SAS tasks ar- 
fgen and rmfgen, respectively. Moreover, source spectra were 
binned to a minimum of 25 counts per bin. 

3. Temporal analysis 

The background subtracted 0.3-10 keV Swift XRT and 
XMM-Newton light curves of the X-ray afterglow of GRB 
050730 are shown in Fig. [2] The same figure also shows the 
BAT 20-150 keV prompt emission light curve, converted in the 
0.3-10 keV energy band using the BAT spectral best fit model 
which is valid also in the XRT bandpass (see Sect. 14. U . The 
UVOT optical light curves in the V and B bands, in arbitrary 
units, are also plotted (see Sect. 12.31 ). 
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Fig. 3. Swift XRT 0.3-10 keV light curve of the GRB 050730 
X-ray afterglow during the first orbit. The solid line is the best 
fit model to the data obtained considering a linear rise exponen- 
tial decay for the three flares (see Sect. 13. 11 1. The dashed vertical 
lines delimit the seven time intervals considered for the spectral 
analysis. Data are binned to 8 seconds resolution and errors are 
at the lcr level. 



3.1. X-ray afterglow 

The X-ray afterglow of GRB 050730 is characterized by 
a very complex structure. The first Swift orbit (from 7o+133 s 
to 7o+794 s), after an initial steep decay phase that joins well 
with the end of the BAT prompt emission, is dominated by three 
bright X-ray flares peaking at about 235, 435 and 685 seconds 
after the BAT trigger. Taking into account cosmological time 
dilation these times correspond to about 47, 88 and 138 sec- 
onds in the GRB rest frame. A flaring episode is also observed 
in the second orbit peaking at about 7o+4500 s. While the un- 
derlying decay of the afterglow during the first two orbits is 
shallow, starting from the third orbit (To +10 ks) the afterglow 
light curve shows a much steeper decline with superimposed 
flaring activity. 

We first modeled the X-ray light curve of the afterglow 
with a double broken power-law model with slopes a\, az, 
and temporal breaks t\ , t%, describing the underlying power-law 
decay of the afterglow, plus seven Gaussian functions model- 
ing the flaring episodes. We found for the first power-law a 
decay index ct\ — 2.1 + 0.3 followed, after a first time break at 
1 1 = 237 + 20 s, by a shallower decay with index — 0.44 
A second temporal break is found at ?2 = 10.1^ ks after 
which a steep decay with index = 2.40*^ is observed. 
This model did not provide a good fit (x 2 r = 1.73, 143 d.o.f.), 
mostly due to short time scale fluctuations and to deviations of 
the three first bright flares from a symmetric Gaussian shape. 

We thus considered a different functional form for the three 
X-ray flares, namely a linear rise exponential decay: F(t) oc 
(f - ?o)/(f p - fo) for times between the flare start time to and 
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Fig. 4. Swift XRT 0.3-1.5 keV (upper panel) and 1.5-10 keV 
(middle panel) light curve of the GRB 050730 X-ray afterglow 
during the first Swift orbit. In the lower panel the corresponding 
hardness ratio is plotted. Data are binned to 12 seconds reso- 
lution and error bars indicate statistical uncertainties at the lcr 
level. 

peak time f p , F(f) oc exp[-(f - t p )/t c ] for t > t p where f c is the 
exponential decay time. The model improves significantly the 
fit with^ = 1.43 (140 d.o.f.) and F-test chance probability of 
1.3 x 10~ 6 . The best fit parameters of the overall underlying 
double broken power-law model were unchanged with respect 
to the previous fit listed above. The linear rise exponential de- 
cay best fit parameters for the first three flares are given in Table 
|2]while Table [3]reports the Gaussian best fits for the other four 
flares. As Fig. [3] illustrates, the second brightest flare (referred 
to as flare 2 in the following) has a flux variation of amplitude 
~ 3.6 and is characterized by a steep rise, lasting ~ 90 seconds, 
followed by a slower decay with duration ~ 170 seconds. An 
asymmetric shape, with a steep rise followed by a shallower 
decay, was also observed for the other two flares (flare 1 and 3) 
and for these episodes the flux variation was ~ 1.7 and ~ 2.6, 
respectively. 

The ratio between the duration (Af) and the peak time (f p ) of 
the X-ray flares was calculated using the linear rise exponential 
decay and Gaussian best fit parameters for the first three and 
the last four flares, respectively. For both analytical models the 
duration of a flare was defined as the time interval during which 
the flare intensity was above 5% of its peak value. The results 
are given in Tables^andlJ] 

The rise and decay portions of the first three flares were 
also fit with single power-laws. For the estimation of power- 
law indices, times were expressed from the onset of the flares 
(?o) for the rise, from the peak time (f p ) for the decay, and the 
contribution of the underlying afterglow decay was taken into 
account. The power-law temporal indices for the rising (a r ) and 
decaying (ad) portions of the three flares are reported in Table 
HI 



6 



M. Perri et al.: The exceptionally extended flaring activity in the X-ray afterglow of GRB 050730 



The temporal behavior of the afterglow during the first three 
X-ray flares was also studied in different energy bands. In 
Fig. |4] we show the 0.3-1.5 keV (upper panel) and 1.5-10 
keV (middle panel) light curves of the early X-ray afterglow 
together with the corresponding hardness ratio (lower panel). 
From the figure it is apparent that i) in the harder band the pro- 
file of the flares is sharper and ii) a peak time shift, with flares 
in the hard band peaking at earlier times, is observed. We note 
that these temporal properties as a function of energy are the 
same observed for the prompt emission pulses of GRBs. 

As already mentioned, late flaring activity of the X-ray af- 
terglow was observed at ~ 4.5, 10.4, 18.7 and 41.2 ks after 
the trigger. Due to the Swift orbital gaps the temporal cover- 
age of flares 4, 5 and 6 is poor and these episodes could not 
be well studied, but the last flare was entirely covered by the 
XMM-Newton follow-up observation allowing a detailed tem- 
poral analysis. In Fig. |5] the XMM-Newton (MOS 1 +MOS2) 
light curve is shown together with the XRT curve covering the 
same time interval. A very good agreement between the two 
curves is found and the ~ 5% higher normalization of the XRT 
data points is of the same order of the uncertainties in the abso- 
lute flux calibration of the instrument (Campana et al. 2006). 

We first fit the XMM-Newton light curve with a single 
power-law decay and obtained an extremely poor fit to the 
light curve, with a fit statistic of x 2 r — 8.81 (32 d.o.f.) and 
a decay index of a = 2.10 + 0.04. The light curve was then 
parameterized with a long duration flare super-imposed on a 
steep power-law decay. We obtained an underlying decay in- 
dex a = 2.45 +0.15, while the flare could be adequately mod- 
eled with a Gaussian function; the fit statistic was acceptable 
(x 2 r = 1.22, 29 d.o.f.). The long duration flare peaked at 41 ks 
(or 8.2 ks in the GRB rest frame), with a cr of 6.8 ks (1.4ks 
rest frame). The total fluence of the flare when parameterized 
this way was 8.3 x 10 50 erg (1.5-50keV band in the GRB rest 
frame), which represents 20% of the integrated afterglow emis- 
sion over the XMM-Newton observation. The rise and decay 
phases of flare 7 were also fit with single power-laws express- 
ing times from the onset and from the peak time of the flare, 
respectively. The measured temporal indices a r and a& are also 
reported in Table [4] 

3.2. Optical band 

The UVOT optical light curve of the afterglow of GRB 
050730 in the V and B bands is also characterized by a complex 
behavior, as is illustrated in Fig. [2] 

The early (r +100 s - r +800 s) UVOT V light curve re- 
veals flaring activity. A re-brightening is in fact observed at 
~ 7o+500 s, almost simultaneously with the brightest X-ray 
flare observed with XRT. Due to the relatively poor sampling of 
the UVOT light curve, a detailed temporal analysis is not possi- 
ble, however we note that the amplitude of the flux variation in 
the V band (~ 3) is of the same order of the one measured in the 
X-ray energy band. The decay of the afterglow in the V band, 
up to about 7o+12 ks and excluding the re-brightening episode, 
is shallow (ffy ~ 0.3) and in agreement with the one measured 
in the X-ray band (ax = 0.44^ ^). Hints of a steepening of 




5x10 4 

Time since BAT trigger (s) 



Fig. 5. XMM-Newton (MOS 1 +MOS2, filled circles) and Swift 
XRT (open circles) 0.3-10 keV light curve of the late (~ 7/ +12 
hours) X-ray afterglow of GRB 050730. The solid line is the 
best fit model to the XMM-Newton light curve. The dashed line 
represent the underlying power-law decay (see Sect. 13. U . Data 
points errors are at the lcr level. 

the V light curve at ~ To+12 ks are also found, but the limited 
temporal coverage and the large statistical uncertainties char- 
acterizing the UVOT data points do not allow us to constraint 
the late optical afterglow decay. 

The UVOT light curve in the B band is characterized by a 
flat power-law decay (a^ ~ 0.3) without re-brightening events. 
However, it should be noted that i) the curve is very sparsely 
sampled, with only three data points during the first Swift orbit, 
and ii) in the B band there is a strong flux reduction due to the 
Lyman break and thus large statistical uncertainties affect the 
data points. 

Accurate optical observations of the afterglow of GRB 
050730 in the R and / bands are presented by Pandey et 
al. (2006). The authors report an early time flux decay with 
indices a R = 0.54 + 0.05 and an = 0.66 ±0.11 followed, after a 
temporal break at about 9 ks after the trigger, by a strong steep- 
ening with decay indices = 1.75 + 0.05 and a\ = 1.66+0.07. 

4. Spectral analysis 

For the spectral analysis of Swift XRT and XMM-Newton data 
we used the XSPEC package (v. 1 1.3.2, IArnaud 19961) included 
in the HEAsoft package (v. 6.0.4). 

4.1. XRT 

As a first step, the 0.3-10 keV XRT average spectrum 
during the first Swift orbit (WT mode, from r +133 s to 
7q+794 s) was fit adopting a single power-law model with the 
neutral hydrogen-equivalent absorption column density fixed 
at the Galactic value in the direction of the GRB (NS = 
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3.0 x 10 20 crrT 2 , |Dickey & Lockman 1990| >. The fit obtained 
was poor (x^ - 1-28, 311 d.o.f.). From the inspection of 
residuals a strong deficit of counts at low energies, likely due 
to the presence of an absorbing column density in excess to 
the Galactic one, was found. Indeed, the addition of a col- 
umn density using Solar metallicity redshifted to the rest 
frame of the GRB host (z = 3.969) free to vary in the spec- 
tral fit (zwabs model in XSPEC) resulted in a more acceptable 
fit with a column density of = (1.28^^) x 10 22 cnT 2 (see 
Table [5]). This value is consistent with the neutral hydrogen 
column densities derived from the optical spectra reported by 
|Starling et al. (2005)| and |Chen et al. (20"05a)l 

From Fig. [4] it is apparent that strong spectral evolution 
takes place during the intense flaring activity observed in the 
first Swift orbit. We thus split the WT observation in seven 
time intervals to study the spectra during the rise and the de- 
cay portions of each flare (see Fig.0. The time-resolved spec- 
tral best fits (see Table [5]l clearly show evidence for spectral 
variation during the flares and an overall softening of the spec- 
tra with time associated with a decrease of the rest frame col- 
umn density. The time-resolved spectral analysis was also per- 
formed adopting a broken power-law model to investigate the 
possible presence of spectral breaks (e.g. Falc one et al. 2 006 ; 
IGuetta ef al. 2007). Also in this case an additional absorption 
column density at the rest frame of the GRB host was con- 
sidered. For all segments we did not find evidence for spectral 
breaks within the XRT energy band. 

The XRT late 0.3-10 keV spectrum (PC mode, from 
7o+4.0 ks to To+ 143.4 ks) was also fit using a single power- 
law model with a fixed Galactic absorption column density and 
an additional absorbing column at the burst rest frame. The ob- 
servation was divided in two time intervals (4.0-18.1 ks and 
2 1.3-143 .4 ks from the trigger). The results are listed in Table[5] 
where we see that spectral softening between the two time in- 
tervals was found. 



4.2. XMM-Newton 

The PN and MOS 0.3-10 keV spectra were fit jointly al- 
lowing the cross normalization between the detectors, which is 
consistent within < 5%, free to vary . The two MOS spectra 
and responses were combined to maximize the signal to noise, 
after first checking that they were consistent with each other. 

As in the case of the XRT spectrum, the PN and MOS spec- 
tra were first fit with a single power-law model with the neu- 
tral hydrogen-equivalent absorption column fixed at the known 
Galactic value. The fit obtained was not acceptable (j^ 2 = 1.40, 
489 d.o.f.), while the energy index was /3 = 0.76 + 0.02. 
The addition of a neutral absorption column in the GRB host 
galaxy frame at z = 3.969 resulted in more acceptable fit 
Or 2 = 1.14, 489 d.o.f.) with an excess absorption column above 
the Galactic column of (6.8 + 1.0) x 10 21 cm 2 , while the con- 
tinuum energy index was now f3 = 0.87 + 0.02. 

The XMM-Newton afterglow spectra were also sliced into 
three segments of approximately lOks in duration, in order to 
search for any spectral evolution within the XMM-Newton ob- 
servation. A small change in the continuum parameters was 



found, the spectrum evolved from hard to soft; the energy 
index changed from f3 = 0.87 + 0.03 in the first lOks, to 
/3 = 0.99 + 0.05 during the final segment. No evidence was 
found for a change in the column density, which was subse- 
quently fixed at A^h = 6.8 x 10 21 cirT 2 in all the segments. The 
spectral best fit parameters are shown in Table[5] 

We also searched for any evidence of emission lines, either 
in the mean spectrum, or in the three segments. No statistically 
significant lines were found, at the level of > 99% confidence. 
As the redshift of the burst is known, then we can set an upper- 
limit to the equivalent width of any emission lines. Over the 
range of 0.4-8 keV and using the mean spectrum, we found a 
< 100 eV upper limit to any emission lines. More specifically 
we can set a limit on the iron Ka line (e.g. the H-like line at 
6.97 keV rest frame, 1.40keV observed frame) of < 30 eV. 

5. Discussion 



5. 1 . Early X-ray light curve 

This GRB, with its rather high redshift, z = 3.969, gives 
us the possibility to investigate the X-ray and optical light 
curves soon after the trigger and thus to study in detail the 
soft tail of the prompt emission and the very beginning of 
the afterglow phase. Indeed, due to cosmological time dilation 
the XRT and UVOT observations started, in the rest frame of 
the burst, only 27 seconds after the trigger. Other examples of 
such early observations of high redshift bursts are GRB 050904 
(z = 6.29,ICusumano et al. 2006bT l and GRB 060206 (z = 4.0, 
IMonfardini et al. 20061 

The XRT light curve shows at the very beginning (133- 
205 s from the BAT trigger) a rapidly decaying emission that 
joins quite nicely with the BAT flux when converted to the 
XRT bandpass, a feature that has been observed in various 



Swift bursts (Tagliaferri et al. 2005 Barthelmy et al. 2005b 
INousek et al. 20061 IQ'Brien et al. 2006t . For GRB 050730. 
the steep temporal decay index {ct\ = 2.1 + 0.3, see 
Sect. I3.lt characterizing the early X-ray light curve sug- 
gests that the observed emission is likely associated with 
the tail of the prompt emission rather than to the beginning 
of the shocked inter-stellar medium afterglow phase. Indeed, 
in the internal shock model scenario for the prompt y-ray 
emission (e.g. IRees & Meszaros 1994t . the cessation of the 
emissivity is characterized by a rapid decay due to the de- 
layed arrival of the high angular latitude prompt emission 
of the shocked surface (high-latitude or curvature emission, 
IKumar & Panaitescu 20001 IDermer 2004i l. The high-latitude 
emission predicts a relationship between the temporal decay in- 
dex a& and the spectral index fid during the decay given by the 
equation — 2+fid (Kumar & Panaitescu 2000) where the de- 
cay slope is measured using as zero time (?o) the beginning of 
the prompt emission (F{t, v) °c (t - to)~ ai v~^ d ). We have tested 
this prediction by fitting the early XRT light curve (segment 
"0" in Fig.0 with the above constraint using the spectral index 
measured during the decay (fj^ = 0.42 + 0.08, see Table|5]l and 



leaving the zero time as a free parameter (see Liang et al. 2006 
for a description of the method). We found to = —31^ s, i.e. 
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the zero time is located at the rising segment of the BAT prompt 
emission light curve (see Fig.[T|), as expected in the framework 
of the high-latitude emission. This result indicates that the steep 
decay observed in the early XRT light curve is most likely the 
soft tail of the BAT emission. 

This hypothesis is also supported by the BAT and early 
(segment "0" in Fig. [3]) XRT spectra. We found very similar 
spectral slopes (y3 B AT = 0.5 ± 0.1 and fi XRT = 0.42 ± 0.08, 
respectively), likely indicating that the early X-ray and y-ray 
emissions are produced by the same emission mechanism. 

5.2. Flaring activity 

The exceptionally extended flaring activity of the X-ray af- 
terglow of GRB 050730 allows us to study this phenomenon 
over more than two orders of magnitude in time. The XRT 
early steep decay is followed by three bright X-ray flares peak- 
ing at 236, 437 and 685 seconds after the BAT trigger (see 
Fig. |3). These flares were under the sensitivity of the BAT 
instrument and thus were not detected in the hard X-ray en- 
ergy band. These flares show, as other strong flares observed 
with Swift (e.g. IRomano et al. 20061 IFalcone et al. 2 006 ; 



Paga ni et al. 2006[IGodet et al. 20061 IChincarini et al. 20061 ), a 



clear spectral evolution with the hardness ratio that mimics 
the variation of the light curve. A phase lag, with the harder 
(E > 1.5 keV) light curve peaking at earlier times with respect 
to the softer energy band is also found (see Fig . |4j . 

So far, most of the X-ray flares observed in Swift X- 
ray light curves have been interpreted as due to late in- 
ternal shocks (e.g. IBurrows et al. 2005bl [Romano et al. 2006; 
Falc one et al. 2 006 ; G odet et al. 20061) . In this scenario the 
GRB central engine is active far beyond the end of the 
prompt y-ray emission phase requiring new mechanisms ca- 
pable of powering new relativistic outflows at late-time (e.g. 
King et al. 2005 IPerna et al. 20061 ). A diagnostic to check if 



the re-brightenings are due to late internal shocks has been re- 
cently proposed by Liang et al. (2006). In the internal-origin 
scenario for X-ray flares the decay emission of the flaring 
episodes should be dominated by high-latitude emission with 
the decay temporal index related to the decay spectral index as 
oid = Pi + 2 dKumar & Panaitescu 20001 ). We have checked this 
hypothesis for the first three bright X-ray flares observed in the 
GRB 050730 afterglow. The spectral indices measured during 
the decays were used (segments "lb", "2b" and "3b", see Fig. [3] 
and Table [5]) and we found best fit zero times fo,i =193*?j s, 
fo,2 =341^2 s and fo,3 =592+„f s, respectively. We can see that 
for all three episodes the zero time values are located at the be- 
ginning of the corresponding flare and thus the observed decay 
slopes are consistent with being due to high-latitude emission 
as predicted by the late internal shock scenario. 

Another possibility is that the X-ray flares are pro- 
duced by refreshed shocks dRees & Meszaros 1 998; 
IKumar & Piran 20001 ISari & Meszaros 20001 ). In the standard 
internal-external fireball model (e.g. Rees & Meszaros 1994; 
ISari & Piran 19971 ) re-brightening episodes in the afterglow 
light curve are explained as slow shells ejected from the 
central engine during the prompt phase that catch up with the 



main afterglow shock after it has decelerated in the external 
inter-stellar medium. Indeed, Guetta et al. (2007) have recently 
interpreted the X-ray flares of the afterglow of GRB 050713A 
as due to refreshed shocks. Here we applied the diagnostic 
proposed for GRB 0507 13A to the case of GRB 050730. First, 
to estimate the rise and decay slopes of the early flares we 
operationally selected as zero time the beginning and the peak 
of the flare, respectively. This assumption implicitly means 
that the flare is completely independent from the main event 
generating the prompt emission and the forward shock light 
curve. We found temporal indices in the 0.8-1.8 range (see 
Table |U). Much steeper slopes (~ 3-6) are obtained if instead 
times are referenced from the BAT trigger, as usually done 
for most of the Swift X-ray flares (e.g. IBurrows et al. 2005bl 
IRomano e t al. 2006; Falcon e~et al. 20061 ) . Second, we checked 
if the temporal and the spectral indices during the decay 
phase are related as predicted by the standard afterglow model 



(Sari, Piran & Narayan 1998). We restricted the analysis to 
the episodes which have enough statistics to allow a detailed 
temporal and spectral study, i.e. the first three flares and the last 
one. For the decay phases we expect the following relations 
to hold ( |Sari Piran & Narayan 19981 |Dai & Cheng 2001| ): 
F v oc (v m lv c T l l l {vlv m Y pl2 F v , nmx cx ( f 0r p > 2 ) 

and F v oc (v m /vJP- r >l 2 (v/v c )-Pl 2 F v , max oc r (? P+ w)n6 (for 
1 < p < 2). From the best fit spectral energy indices during 
the decay of the four flares (see Table [5]) we derived the 
predicted temporal slopes a\ = 0.93 + 0.05, a\ = 0.89 + 0.03, 
a\ = 1.02 + 0.15, and a\ = 1.40 + 0.05, which are consistent 
with the measured values of the decay indices listed in Table 
|U We note that for the last flare, peaking at 41.2 ks after the 
trigger, we used the afterglow model relations for the slow 
cooling case (F v oc {v/v m y^ l ^ 2 F yjnax oc fX 1- ?)/ 4 ) which very 
likely applies at these late times (Sari, Piran & Narayan 1998). 

Important information on the mechanism producing the 
flares can also be obtained from the comparison between the 
observed variability timescale and the time at which the flare 
is observed dloka et al. 2 005). We calculated, for all the 7 re- 
brightenings of the GRB 050730 afterglow, the ratio between 
the duration of the flare (Af) and the time when the flare peaks 
(f p ) (see Tables [2] and EJ. We found At/t p ~ 0.3 for the first 
flare and Af/f p ~ 0.6 - 0.7 for the others, in agreement with 
the Af/fp > 0.25 limit discussed by Ioka et al. (2005) for the 
refreshed shock scenario. Remarkably, the Af/f p ratio is nearly 
constant for all flares, with the exception of the first one which 
occurred during the bright and steep tail of the prompt emis- 
sion (see Fig. [3} and most likely we are observing only the tip 
of the flare. A duration of the flare proportional to f p , as ob- 
served in the afterglow of GRB 050730, is explained both in 
the refreshed shock scenario (IKumar & Piran 20 00 ) and in the 
late internal shock model ( IPerna et al. 20061 ). We note that the 
last four flares have Af/f p values in the range 0.7-0.9, while 
for flares 2 and 3 we observe slightly lower ratios (0.5-0.6), 
possibly indicating a moderate temporal evolution of Af/f p . 
However, due to the relatively large error bars, the measured 
values are also consistent with being constant and ~ 0.6 - 0.7. 

Flaring activity is also observed in the optical afterglow 
of GRB 050730. Indeed, the UVOT V early light curve (see 
Fig. [2]) shows a re-brightening at ~ 7q+500 s, almost simul- 
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taneously with the brightest X-ray flare observed with the 
XRT peaking at 7o+437 s. Moreover, the optical re-brightening 
amplitude (~ 3) is of the same order of the flux variation 
observed in the X-ray energy band. Strong indications of 
correlated variability between the X-ray and optical energy 
bands are also present in the NIR/Optical light curves re- 
ported by Pandey et al. (2006): a significant re-brightening 
at about 7o+4 ks (R filter), close to the X-ray flare peaking 
at To+4.5 ks (see Fig. |5J, is observed. Moreover, a bump in 
the J, I, R and V light curves at about 7o+10 ks is observed 
( Pandey et al. 2006} , again simultaneous with the X-ray flare 
peaking at To+10.4 ks. Although the X-ray flaring activity is 
not uniformly covered by optical observations, we find several 
indications of simultaneous re-brightening events in the X-ray 
and optical bands, in agreement with the refreshed shock model 
(e.g. lGranot etal. 2003b . 



5.3. Evidence of a jet break 

The X-ray afterglow light curve shows a clear temporal 
break around 10 ks after the trigger (see Sect. 13. lb . At about 
the same time a strong steepening of the / and R afterglow 
light curves, although with a shallower post-break slope, is also 



observed (see Pandey et al. 2006 for a detailed analysis on the 
discrepancy between the X-ray and optical slopes). Due to the 
achromatic nature of the temporal break it is very likely that a 
jet break is occurring when the bulk Lorentz factor y of the col- 
limated relativistic outflow becomes lower than the inverse of 
the jet opening angle 0j et (e.g., IRhoads 19971 ISari et al~9 99 ), 
as also reported by Pandey et al. (2006). In this framework, 
the jet opening angle can be determined through the equation 
jet = 0.161[f b /(l + z)] 3/8 («/7/£ is o) 1/8 (e.g.. Bloom etal. 20031) 
where 9 jel is in radians, the jet temporal break fb in days, the 
total isotropic-equivalent energy E- lso in units 10 52 erg, the den- 
sity n of the circumburst medium in ctrT 3 and r\ is the efficiency 
of conversion of the outflow kinetic energy in electromagnetic 
radiation. 

An accurate estimation of the bolometric isotropic- 
equivalent energy radiated by GRBs requires the knowl- 
edge of their intrinsic spectrum over a broad energy band 
dBloom et al. 20011 lAmati et al. 20021 . For GRB 050730, the 
spectrum is well fit by a single power-law with energy index 
Pbat = 0.5 +0.1 up to the high energy limit of the BAT sen- 
sitivity bandpass (150 keV) indicating that i) we are observing 
the low energy tail of the Band model (B and et al. 19 93 ) gen- 
erally used to describe GRB spectra and ii) the vF(v) spectrum 
peak energy E p is above ~ 750 keV (i.e. (1 + z) x 150 keV) in 
the burst rest frame. A lower limit to the bolometric isotropic- 
equivalent radiated energy Ei so is given by the observed radi- 
ated energy in the BAT bandpass (E^ T = (8.0+ 1.0) x 10 52 erg, 
see Sect. 12. U . An upper limit to E[ so is obtained in the most 
conservative case where the peak energy E p is equal to 10 4 keV 
( Am ati et al. 2002 ). By integrating the best fit BAT power-law 
spectrum in the whole 1-10 4 keV rest frame energy band we 
thus derived an upper limit of 4.5 x 10 53 erg to Ei so . Taking 
the central value of the interval derived above we obtained for 
GRB 050730 £ iso = (2.6 ± 1.9) x 10 53 erg. 



With z = 3.969, % = 10. ks (see Sect. ED and the E i$0 
range above derived, we find #j et = 1.6^2*| ^ e §' f° r H = 0.2 
dFrail etal. 20011) and assuming the value of circumburst den- 
sity n — 10 ctrT 3 discussed by Bloom et al. (2003). With 
this value of the jet opening angle, the inferred collimation- 
corrected bolometric radiated energy is £j et = (LO^ng) x 
10 50 erg. 

Taking into account the E\ so and £j et values derived above 
and the lower limit to the peak energy (E v > 750 keV rest 
frame), we find that GRB 050730 is consistent with the E v 
vs. E^ relation found by Amati et al. (2002) and recently 
updated in Amati (2006). We also find that GRB 050730 is 
inconsistent, even taking into account the 3<x scatter around 
the best fit correlation, with the E p vs. £j et relation found by 
Ghirlanda et al. (2004) and with its updated version presented 
by Nava et al. (2006). In order to make this GRB consistent 
with the £p-£j et relation a much higher circumburst density 
(n ~ 10 5 cm 4 ) would be required. GRB 050730 is inconsis- 
tent with the model-independent E^-Ep-tb correlation found 
by Liang & Zhang (2005). 

6. Summary and conclusions 

We have presented a detailed temporal and spectral analy- 
sis of the afterglow of GRB 050730 observed with Swift and 
XMM-Newton. The most striking feature of this GRB is the in- 
tense and exceptionally extended, over more than two order of 
magnitude in time, X-ray flaring activity. 

Superimposed to the afterglow decay we observed seven 
distinct re-brightening events peaking at 236 s, 437 s, 685 s, 
4.5 ks, 10.4 ks, 18.7 ks and 41.2 ks after the BAT trigger. The 
underlying decline of the afterglow was well described with a 
double broken power-law model with breaks at t\ = 237 + 20 s 
and ?2 = 10.1^2 2 ks. The temporal decay slopes before, be- 
tween and after these breaks were cr\ =2.1 +0.3, cti - 0.44^j ( 'g 
and q-3 = 2.40^qq^, respectively. 

Strong spectral evolution during the flares was present to- 
gether with an overall softening of the underlying afterglow 
with the energy index varying from (3 - 0.42 + 0.08 dur- 
ing the early (133-205 s) steep decay to ft = 0.99 + 0.05 
at much later (50-60 ks) times. An absorbing column density 
/Yy = (1.28^- 2 g) x 10 22 cnT 2 in the host galaxy is observed 
during the early (133-781 s) Swift observations while a lower 
column density N^. = (0.68 + 0.10) x 10 22 cm~ 2 is measured 
during the late (29.4-50.8 ks) XMM-Newton follow-up obser- 
vation, likely indicating photo-ionization of the surrounding 
medium. Evidence of flaring activity in the early UVOT op- 
tical afterglow, simultaneous with that observed in the X-ray 
band, was found. 

From the temporal analysis of the first three bright X-ray 
flares we found that the rise and decay power-law slopes are 
in the range 0.8-1.8 if the beginning and the peak of the flares 
are used as zero time, respectively. We also found that, with the 
exception of the first flare, for all episodes the ratio between the 
duration of the flare (At) and the time when the flare peaks (f p ) 
is nearly constant and is Af/f p ~ 0.6 - 0.7. 

We showed that the observed properties of the first three 
flares are consistent with being due to both high-latitude emis- 



10 



M. Perri et al.: The exceptionally extended flaring activity in the X-ray afterglow of GRB 050730 



sion, as expected for flares produced by late internal shocks, 
or to late time energy injection into the main afterglow shock 
by slow moving shells (refreshed shocks). An analysis of a 
larger sample of bursts would help in understanding what are 
the physical mechanisms responsible for the X-ray flaring ac- 
tivity. 

We interpreted the X-ray temporal break at around 10 ks 
as a jet break and derived a cone angle of ~ 2 deg and a radi- 
ated energy £j et = (0.2 - 3.3) x 10 50 erg against an isotropic- 
equivalent energy E iso = (0.7 - 4.5) x 10 53 erg. GRB 050730 
satisfies the E p vs. E lso Amati relation while is inconsistent with 
the E p vs. £jet Ghirlanda relation. 
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Table 1. UVOT detections of GRB 050730 in the V and B filters. Column (1) gives the image mid time in seconds since the BAT 
trigger, column (2) the net exposure time, column (3) the filter used, and column (4) the afterglow magnitude with lcr error. 



T(mid) 


Exposure 


Filter 


Magnitude 


(s) 


(s) 






170 


99 


V 


17.4 ±0.2 


270 


10 


B 


18.4 + 0.4 


299 


9.7 


V 


17.2 ±0.5 


383 


9.7 


V 


17.3 ±0.5 


468 


9.7 


V 


16.4 ±0.3 


524 


30 


B 


19.0 ±0.4 


552 


9.7 


V 


16.3 ±0.3 


637 


9.7 


V 


17.0 ±0.4 


721 


9.7 


V 


16.8 ±0.4 


734 


20 


B 


19.2 ±0.5 


10164 


900 


B 


20.2 ± 0.3 


11947 


837 


V 


18.6 ±0.2 


23519 


835 


V 


19.5 ±0.3 


34990 


843 


V 


19.9 ±0.4 



Table 2. Temporal best fit parameters of the first three bright X-ray flares of GRB 050730 using a linear rise exponential decay 
model. The corresponding Af/f p values for the three flares are also indicated (see Sect. 13. 11 1. 



Flare 


to 


h 


fc 


K x lO" 1 " 


A'/fp 




(s) 


(s) 


(s) 


( erg cirr 2 s -1 ) 




1 


207^ 


236+j 


n +b 5 


8 + 1 


0.3 ±0.1 


2 
3 


344+5 
6141' 


437!= 
685!? 


54^ 


12.2 ±0.7 
6.1 ±0.7 


0.58 + 0.04 
0.5 ±0.1 



Table 3. Temporal best-fit parameters of the X-ray flares 4, 5, 6 and 7 of GRB 050730 using Gaussian functions. An asterisk 
indicates a frozen parameter. The corresponding Af/f p values are also indicated (see Sect. 13. U . Best-fit parameters for flare 7 
were derived from XMM-Newton data only. 



Flare 


tp 


a 


a:x io -11 


Af/'p 




(s) 


(s) 


( erg cnr 2 s~') 




4 

5 


44841™ 
10391 + ^o 


641+J59 
wl -149 

1500* 


11 ±2 
5.5 ± 1.0 


0.7 ± 0.2 
0.71 ±0.02 


6 
7 


-744 


3638 + ^gg 
6170!- 


1.5 ±0.4 
0.12 ±0.02 


0.9 ±0.2 
0.7 ±0.1 



Table 4. Best fit temporal indices of the rising (a r ) and decaying (ofj) portions of the X-ray flares 1, 2, 3 and 7 using a single 
power-law model. 



Flare 


a c 






1 


-1.56 ±0.69 


1.25 


±0.32 


2 


-1.84 ±0.30 


1.31 


±0.44 


3 


-1.46 ±0.62 


0.75 


±0.35 


7 


-0.89 ± 0.34 


1.61 


±0.44 
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Table 5. Results of single power-law spectral fits to the 0.3-10 keV spectrum of the afterglow of GRB 050730. A local (z 
absorption column fixed at the known Galactic value of iV§ = 3.0 x 10 20 cirT 2 (Dickey & Lockman 1990) was used in the 
An asterisk indicates a frozen parameter. 



segment 


time interval 

(s) 


N~ H x 10 22 
(cm" 2 ) 


j8 


Xr (d-O.f.) 


WT (all) 


133-781 


1.28«g 


70 +0 03 

' -0.03 


1.01 (310) 



WT(0) 


133-205 


1.8 


WT (la) 


205-233 


1.6 


WT (lb) 


233-313 


3.1 


WT (2a) 


313-433 


2.1 


WT (2b) 


433-601 


0.9 


WT (3a) 


601-681 


0.7 


WT (3b) 


681-781 


1.0 


PC(1) 


4001-18149 


1.1 


PC (2) 


21288-143438 


1.0 



+0.9 
-0.8 
+2.5 
-1.6 
+1.3 
-1.1 
+0.8 
-0.7 
+0.5 
-0.5 
+0.S 
-0.7 
+0.6 
-0.6 



+0.4 
-0.4 
+0.7 
-0.6 



47+0.08 

"• HZ "-0.08 





86 


(76) 







86 


(24) 


0.82 + «|? 


1 


30 


(51) 


07 1+0.08 
"•' 1 -0.08 


1 


10 


(93) 


70+ 07 

u - /u -0.07 


1 


11 


(111) 


77+°- 12 

-0.12 





87 


(46) 


i or 010 

1 - ul -o.io 





81 


(65) 



61+ 04 

"• ul -0.04 

n ft i + 08 

u - 8 1-0.08 



0.95 (224) 
1.16 (86) 



XMM(all) 29436-59811 0.68! 



0.87!»;JJ 2 



1.14 (489) 



XMM (1) 


29436^10000 





68* 





O-7+0.03 
8 '-0.03 


0.88 


(345) 


XMM (2) 


40000-50000 





68' 





q i+o.03 

7 -0.03 


0.98 


(243) 


XMM (3) 


50000-59811 





68* 





qq+o.05 

7 -0.05 


0.77 


(135) 



